Polarization of prompt J /if) in proton-proton collisions at RHIC 



Hee Sok Chung and Chaehyun Yu 
Department of Physics, Korea University, Seoul 136-701, Korea 

Seyong Kim 

Department of Physics, Sejong University, Seoul 143-747, Korea and 
School of Physics, Korea Institute of Advanced Study, Seoul 130-722, Korea 



o 

(N 
C 

in 



in: 

CD 
> 



(N 



<3\ 
O 



X 



Jungil Lee 

Department of Physics, Korea University, Seoul 136-701, Korea and 
Korea Institute of Science and Technology Information, Daejeon, 305-806, Korea 

(Dated: January 15, 2010) 

Within the framework of the nonrelativistic QCD (NRQCD) factorization approach, we compute 
the polarization of prompt J/tp produced at the Brookhaven's Relativistic Heavy- Ion Collider from 
proton-proton collisions at the center-of-momentum energy yfs = 200 GeV. The perturbative contri- 
butions are computed at leading order in the strong coupling constant. The prediction reveals that 
the color-singlet contribution severely underestimates the PHENIX preliminary data for the differ- 
ential cross section integrated over the rapidity range \y\ < 0.35 and its contribution is strongly 
transversely polarized, which disagrees with the PHENIX preliminary data. After including the 
color-octet contributions, we find that the NRQCD predictions for both the cross section and po- 
larization over the transverse-momentum range 1.5 GeV < pr < 5 GeV (1.5 GeV < pr < 2 GeV) 
integrated over the rapidity range \y\ < 0.35 (1.2 < \y\ < 2.2) agree with the data within errors. 

PACS numbers: 13.88.+e,13.85.Ni,13.85.-t 



In the mid-1990s, the Collider Detector at Fermi- 
lab (CDF) Collaboration at the Fermilab Tevatron ob- 
served remarkable surplus of the prompt J/tp and tp(2S) 
production at large transverse momenta (pr) [H 0. 
Based on the nonrelativistic quantum chromodynamics 
(NRQCD) factorization approach jH, Braaten and Flem- 
ing explained the large cross section by introducing the 
color-octet gluon fragmentation mechanism, which is the 
NRQCD extension of the color-singlet version proposed 
in Ref. Q. According to the proposal the prompt 
spin-triplet S-wave charmonia are dominantly produced 
from the color-octet spin-triplet S- wave charm-anticharm 
pairs [cc8( 3 <Si)] that are fragmented from the gluons at 
large px ■ 

Using the approximate heavy-quark spin symmetry of 
the NRQCD Lagrangian, Cho and Wise predicted that 
the spin-triplet 5-wave charmonia H must be trans- 
versely polarized at large pt [yj , if the color-octet gluon 
fragmentation is the dominant source in that region. 
Subsequently, phenomenological polarization predictions 
for the prompt tp(2S) 0-Q and J/tp [13, EH at the Teva- 
tron were carried out. The result is that these states are 
the more transversely polarized the larger pt is @~GjI 
and the rise of the transverse polarization for the J / ip is 
slower than that of the tp(2S) because of the feed-downs 
from the P-wave spin-triplet states \ c j for J = 0, 1, and 
2 [13, [HI]- Therefore, the observation of transverse H at 
large pt is an important independent confirmation of the 
color-octet gluon fragmentation mechanism. 

However, the run-I CDF measurement with an inte- 
grated luminosity of about 110pb _1 [l2j was in dra- 
matic contradiction to these predictions. To make mat- 



ters even more confusing, the run-II CDF data of about 
800 pb -1 13] neither agree with the predictions nor their 
own run-I measurement. The prompt J/tp's from run-II 
CDF data are almost unpolarized and even show a mod- 
erate increase in the longitudinal fraction from pt — 5 to 
30 GeV 



13] 



Further theoretical studies include analyses at hi ghe r 
orders in the strong coupling a* for the color-singlet [14| — 
[IH and for the color-octet [l9[ channels and the rela- 
tivistic corrections [13] ■ Even after those corrections, the 
color-octet mechanism still contributes to the production 
rate significantly. Recent color-singlet-model calculations 
of the polarization of direct J/tp at next-to-leading or- 
der (NLO) in a s predict strong longitudinal polariza- 
tion [T3, QjJ , which disagree with the run-II CDF data 
over the whole range of pr- These disagreements have 
cast doubt on our current understanding of the char- 
monium production mechanism and lead one to analyze 
hadroproduction of charmonia at various circumstances. 

Recently, the PHENIX Collaboration at the 
Brookhaven's Relativistic Heavy-Ion Collider (RHIC) 
reported the cross section [2l| and the polarization [22[ 
of the inclusive J/tp production in pp collisions at the 
center-of-momentum (CM) energy y/s = 200 GeV. The 
data also contain the J ftp's from the B decay, which 
occupy a tiny fraction (0.0361q'o2;s) °f the total rate [23| . 
Therefore, the inclusive charmonium cross section should 
be essentially the same as the prompt one. There are 
several theoretical analyses on the J/tp production in 
pp collisions at RHIC. Several years ago, Cooper, Liu, 
and Nayak reported the NRQCD prediction for the 
production rate for prompt J/tp [24], [2f|, which agrees 
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with the run-3 PHENIX data [26|, |27f. The authors 
of Ref. [HI explained the J/tp production rate with 
two phenomenological parameters of a new production 
mechanism so-called s- channel cut [29]. However, the 
s-channel-cut prediction for the polarization disagrees 
with the PHENIX preliminary data for the large-rapidity 
region 1.2 < \y\ < 2.2 while that fo r \y \ < 0.35 is con- 
sistent with the measurement [22l . |3C| . Furthermore, 
Artoisenet and Braaten [3l[ showed that the s-channel 
cut can be identified with the charm-pair rescattering 
and is not a dominant mechanism for charmonium 
production in high-energy collisions. Very recently, 
Brodsky and Lansberg carried out calculations for 
direct J/tp in the color-singlet model at NLO in a s [32T ]. 
including contributions from eg fusion. They obtained 
a rapidity distribution in good agreement with the 
PHENIX data 32] under the assumption that about 
40% of the J/tp events come from higher resonances. 

In this paper, we present a quantitative analysis of 
the cross section and the polarization of prompt J/tp 
produced in pp collisions at y/s — 200 GeV using the 
NRQCD factorization formalism [3j. The perturbative 
contributions are computed at leading order (LO) in a s . 
By including both the color-singlet and color-octet con- 
tributions, we find that the NRQCD predictions for both 
the cross section and polarization agree with the data 
within errors. 

The schematic form of the NRQCD factorization for- 
mula for the differential cross section of the polarized 
charmonium H\ with momentum P and spin quantum 
number A is given by 



(1) 



where the summation is assumed over the index n for 
the color and angular momentum states of the cc pair. 
The short-distance coefficients da CCn ( p \ which are insen- 
sitive to the long-distance nature of the quarkonium, are 
calculable using perturbative QCD. The nonperturbativc 
nature of the H is factorized into the NRQCD matrix el- 
ements {On )■ The matrix elements are, in general, 
tensors depending on P and the polarization tensor of the 
H\. After using the symmetries of NRQCD, one can re- 
duce these polarized matrix elements (On } in terms 
of the scalar matrix elements (Of) that are independent 
of P and A. The numerical importance of the NRQCD 
matrix elements can be estimated based on the velocity- 
scaling rules of NRQCD [3(. For the spin-triplet S'-wave 
quarkonia H — J/tp or tp(2S), the most important ma- 
trix elements are (Of ( 3 <Si)) for the color-singlet state, 
and (Oi^So)), (Of ( 3 Si)), and (Of ( 3 P )) for the color- 
octet states, respectively. Because prompt J/tp's include 
the samples that come from the decay of tp(2S) and Xcj, 
we also have to consider the production of XcJ- For the 
Xcj, the matrix elements (O^°( 3 P )) and (Of°( 3 5i)) 
are equally important. 

In the PHENIX experiment of pp collisions, they probe 
a moderate pt range, where the fragmentation effect is 
negligible 25j. Therefore, the production of H should be 



dominated by the fusion contributions. Then the differ- 
ential cross section for H\ is expressed as 

do 11 ^) = f. /p ® f j/p ® da^ P) (Of ^>>, (2) 

where fi/ p (x,p) is the parton distribution function 
(PDF) and the sums over the partons i and j are implied. 
Here, x and /i are the longitudinal momentum fraction 
of the parton and the factorization scale, respectively. 

We proceed to compute the differential cross section 
for the H\. In order to make a prediction that is con- 
sistent with that for the Tevatron [lj|, we follow the 
strategies of Ref. except that we neglect the frag- 
mentation. We include the parton processes ij — > cc+k, 
with i, j — g, q, q and q — u, d, s, and neglect heavy 
partons like c and b. Corresponding LO parton cross 
sections da of order a 3 are given in Refs. 0, [ill, [HI- 
The numerical values for the relevant NRQCD matrix 
elements are given in Table I of Ref. [lfj. 1 For the 
PDF's, we choose MRST98LO [H as the default value 
and CTEQ5L 35] for comparison. We use the trans- 
verse mass tut — (4m^ + pip) 1 !" 1 as a common scale 
for the factorization scale and the renormalization scale 
with the charm-quark mass m c = 1.5 GeV. We evaluate 
a s from the one-loop formula using the value of Aqcd 
given in each PDF set [§3, l35| . We estimate theoreti- 
cal uncertainties in our numerical calculations following 
Ref. [Ioj |. The errors are from the matrix elements in 
Table I of Ref. [l(| and from the variations of /i and m c 
within the ranges ^m,T~2mT and 1.45-1.55 GeV with 
the central values and 1.5 GeV, respectively. The 
errors in PDF are estimated by taking the difference be- 
tween MRST98LO and CTEQ5L. We put the color-octet 
1 5o and 3 P matrix elements as (Of ( 1 S'o)) = xM^ , 
(Of ( 3 Po))Mc = (1 - x)Mf jr and vary x from to 
1 with the central value \ . All of the errors listed above 
are added in quadrature. 

One of the most convenient measures of the polarized 
cross section for the spin-triplet S'-wave quarkonium H 
is the variable a defined by 



= (a T - 2cr L )/{a T + 2a L ), 



(3) 



where <jt (cl) is the cross section for the transversely 
(longitudinally) polarized H . For the complete trans- 
verse (longitudinal) case, a = +1 (—1). If the H is unpo- 
larized, then a = 0. 2 In the recent PHENIX analysis |22j . 
they used the hadron CM frame like the CDF analy- 
ses [l3, 21]. Employing the method given in Refs. US, 
one can compute the polarized cross sections ot and <tl 



1 The color-octet matrix elements (Of ( 1 S'o)} and (Of ( 3 Pj)) are 
determined not separately but as a linear combination M^ 1 = 
(Of ( 1 So)} + r{Q§ ( 3 Pp))/m^ which depends on the variable r 
given in Ref. [10(| . 

2 The polarized cross section varies depending on the spin quanti- 
zation axis 36. 1371 . 
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FIG. 1: The differential cross section of prompt J/ip in pp 
collisions at ^/s = 200 GeV as a function of pr in units of 
nb/GeV 2 against the PHENIX preliminary data [II] for the 
inclusive J/ip. The shaded band represents the NRQCD pre- 
diction and the band surrounded by a solid curve is the color- 
singlet contribution at LO in a 3 . 



FIG. 2: The polarization parameter a for the prompt J/ip 
in pp collision at ^fs = 200 GeV as a function of pr against 
the PHENIX preliminary data 22] for the inclusive J/ip. The 
shaded band represents the NRQCD prediction and the band 
surrounded by a solid curve is the color-singlet contribution 
at LO in a s . 



in the hadron CM frame from Eq. ([2]). Substituting these 
results to Eq. ([3]), we obtain a. 

Our prediction for the differential cross section for the 
prompt J/ip is shown in Fig.[T]as a function of px against 
the PHENIX preliminary data (22j. 3 The rate is averaged 
over the midrapidity region and its explicit normalization 
is given by 

^(d 2 -/dyd PT ) y ^^^J_ Y ^-dy, (4) 

where Y — 0.35 is the rapidity cut and B ee is the branch- 
ing fraction for J/ip — > e + e~ . The shaded band indicates 
the NRQCD prediction at order a 3 . The dominant source 
of the theoretical uncertainties is from the scale /x, which 
produces errors of about 78% (94%) at px = 1.5 GeV 
(5 GeV). The variation of m c produces errors of about 
11 % (2 %) at p T = 1.5 GeV (5 GeV). The remaining con- 
tributions from the matrix elements, x, and PDF's pro- 
duce errors of about 11 % (4 %) at p T = 1.5 GeV (5 GeV). 



In the region px > 1.5 GeV, the NRQCD predic- 
tion agrees with the data within theoretical uncertain- 
ties. The curve for the central value tends to overesti- 
mate (underestimate) the data in the region pt < 3 GeV 
(> 3 GeV). In Fig. [TJ the color-singlet contribution at 
LO in a s is displayed as a band surrounded by a solid 
curve, which severely underestimates the data over the 
whole range of px- In the lower px region, the fixed-order 
calculation of order a 3 fails to give reliable predictions. 
In order to extend the prediction to the lower px region, 
one must include the order-a^ 2 — > 1 parton processes, 
its NLO contribution of order a 3 , and multiple soft-gluon 
emissions that should be resummed, 4 which are out of the 
scope of this work. 

Our results for the a of the prompt J/ip integrated 
over the midrapidity region \y\ < 0.35 is shown in Fig. [5] 
against the PHENIX preliminary data [22j. The shaded 
band represents the NRQCD prediction and the band 
surrounded by a solid curve is that for the color-singlet 
contribution. The uncertainties of a are computed in the 
same manner as those of the differential cross section in 



3 The PHENIX preliminary data quoted in Fig. [T]do not include 
an additional global systematic error of 10%. 4 See, for example, Ref. [38t |39H and references therein. 



4 



Fig. [TJ In this case, a large portion of the uncertainties 
cancel because a is computed as the ratio of polarized 
cross sections. In contrast to the cross section, the dom- 
inant sources of the errors are the matrix elements, x, 
and PDF's that account for 99 % in total. 

We predict that the prompt J/ip is almost unpolar- 
ized or slightly longitudinally polarized in the range 
1.5 GeV < p T < 5GeV. Although the NRQCD predic- 
tion has large theoretical uncertainties, it agrees with the 
two data points within errors. As we have expected, the 
fragmentation dominance does not occur and, therefore, 
a strong transverse polarization is not observed. The 
color-singlet contribution is strongly transversely polar- 
ized with small uncertainties, which is disfavored by the 
data. 

The polarization variable a for the large-rapidity re- 
gion 1.2 < \y\ < 2.2 was also measured by the PHENIX 
Collaboration at p T = 1.6 GeV as a = 0.02 ± 0.16 [Hj]. 
Our predictions of the color-singlet model and NRQCD 
are a C oior-sin g ict = 0.19±0.03 and a N RQCD = 0.15±0.10, 
respectively. On the other hand, one can compare the 
differential cross section in the large-rapidity range near 
p T = 1.6 GeV. We hnd that /(1.6 GeV) co ior- S ingiet = 
0.08 ± 0.05 nb/GeV 2 and /(1.6 GeV) N R Q CD = 2.24 ± 
1.82nb/GeV 2 , where f(p T ) is the same as the right side 
of Eq. (Q| except that the average is over the range 
1.2 < \y\ < 2.2. The color-singlet model prediction 
severely underestimates the differential cross section and 
the NRQCD prediction agrees with data in both a and 
f(pr) in the large-rapidity region. 

Our results in Fig. [2] are the first NRQCD predictions 
that can be compared with the PHENIX results for the 
Pt distributions of the polarization [22j. Our px distri- 
bution for the production rate shown in Fig. [T] integrated 
over the midrapidity range agrees with a previous result 
in Ref. [25| within errors. Very recently, Brodsky and 
Lansberg computed the rapidity distribution of J/ip in- 
tegrated over pt within the color-singlet model at NLO 
in a s . Unfortunately, we are unable to compare our re- 
sults directly with those in Ref. [32[ because they did not 
calculate the pt distribution and because our LO calcu- 
lation breaks down at small pr so we cannot calculate the 
rapidity distribution integrated over the whole pt range. 



In conclusion, we have provided the NRQCD predic- 
tions for the pt distribution of the differential cross sec- 
tion and the polarization variable a for the prompt J ftp 
produced in pp collisions at y/s = 200 GeV. The short- 
distance processes were computed at the fixed-order a^.. 
We have chosen the numerical values for the nonper- 
turbative NRQCD matrix elements that were fit to the 
prompt J/ip cross section measured at the Tevatron and 
were used to compute the NRQCD prediction for the a 
at the Tevatron. The prediction reveals that the color- 
singlet contribution severely underestimates the data for 
the differential cross section. The NRQCD predictions 
for both the cross section and polarization agree with the 
data within errors over the range 1.5 GeV < pt < 5 GeV 
(1.5 GeV < pt < 2 GeV) integrated over the rapidity 
range \y\ < 0.35 (1.2 < \y\ < 2.2). Our fixed-order cal- 
culation breaks down for pt below 1.5 GeV. In order to 
extend the NRQCD prediction to the lower pt region, 
one must include the 2 — > 1 processes at NLO accuracies 
and resum large logarithms from soft-gluon emissions. It 
would be exciting to see whether the future NRQCD pre- 
diction at this region agrees with the measured PHENIX 
preliminary data at the lowest pt bin. 
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